Dealing with the regional land surface heat fluxes over heterogeneous landscape of desertification area near sea is very important and not an easy problem. In this study,a parameterization method based on Landsat-7 ETM+data and field observations is described and tested for deriving the regional land surface variables, vegetation variables and land surface heat fluxes over heterogeneous landscape.As case study, the method is applied to Jiddah area of Saudi Arabia. The regional distributions of surface reflectance,NDVI, MSAVI,vegetation coverage,LAI,surface temperature, net radiation, soil heat flux, sensible heat flux and  latent heat flux have been determined over Jiddah area. The derived results have been validated by using the“ground truth”. The
resultsshow that the more reasonable regional distributions of land surface variables (surface reflectance, surface temperature), vegetation variables (MSAVIand vegetation coverage), net radiation, soil heat flux and sensible heat flux can be obtained by using the method proposed in this study. Further improvement of the method is also discussed.